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ABSTRACT 

Can we determine a spin parameter of a black hole by observation of a black hole shadow in an accretion 
disk? In order to answer this question, we make a qualitative analysis and a quantitative analysis of a shape and 
a position of a black hole shadow casted by a rotating black hole on an optically thick accretion disk and its 
dependence on an angular momentum of a black hole. We consider two types of inner edges, an event horizon and 
a marginally stable orbit, of accretion disks. We have found black hole shadows with a quite similar size and a 
shape for largely different black hole spin parameters and a same black hole mass. Thus, it is practically difficult 
to determine a spin parameter of a black hole from a size and a shape of a black hole shadow in an accretion disk. 
On the other hand, a displacement of a black hole shadow from a rotation axis of a black hole is caused by a frame 
dragging effect in the vicinity of a black hole. The extent of the displacement largely depend on a black hole spin 
parameter. However, it is difficult to determine a position of a rotation axis of a black hole observationally. So, we 
can not practically measure a spin parameter of a black hole from the displacement from a rotation axis of a black 
hole shadow. We newly introduce a bisector axis of a black hole shadow named a shadow axis. We define this 
shadow axis as a bisector perpendicular to a line segment of the maximum width of a black hole shadow. We can 
determine a position of a shadow axis of by observation of a black hole shadows. For a non-rotating black hole, 
the minimum interval between a mass center of a black hole and a shadow axis is null. On the other hand, for a 
rotating black hole a shape and a position of a black hole shadow are not symmetric with respect to a rotation axis 
of a black hole shadow. So, in this case the minimum interval between a mass center of a black hole and a shadow 
axis is finite. An extent of this minimum interval is roughly proportional to a spin parameter of a black hole for a 
fixed inclination angle between a rotation axis of a black hole and a direction of an observer. The maximum extent 
of these minimum intervals is about \ .5rg. Here, rg is a gravitational radius. This is realized when the maximum 
inclination angle and the maximally rotating black hole in an accretion disk with its inner edge of an event horizon. 
In order to measure a spin parameter of a black hole, if a shadow axis is determined observationally, it is crucially 
important to determine a position of a mass center of a black hole in a region of a black hole shadow. We also 
discuss how to determine a mass center of a black hole by observation of a black hole shadow in an accretion disk. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: nuclei 



1. INTRODUCTION 

During the last century, masses of supermassive black holes 
and stellar mass black holes are observationally measured. De- 
termination of a spin parameter of a black hole is one of the 
greatest challenges in astrophysics in this century. So far, sev- 
eral determination methods of a black hole spin are proposed. 
Broad and/or skewed X-ray Fe Ka emission lines, first discov- 
ered by ASCA, are extensively studied (e.g. Fabian et al. 1995, 
Tanaka et al. 1995, Bromley, Miller & Paiiev 1998). If appro- 
priate accretion disk models in Kerr spacetime are assumed, one 
can explain these broadening features for some spin parameters 
of central black holes. Some authors actually propose some val- 
ues of spin parameters of supermassive black holes (e.g. Miller 
et al. 2002). Melia et al. (2001) said that they constrain the 
black hole spin in Sagittarius A* by sub-millimeter timing ob- 
servations. Recently, by high-resolution infrared observations 
of Sagittarius A* Genzel et al. (2003) estimates the spin param- 
eter of the massive black hole in Sgr A*. These two methods 
crucially depend on a period of a matter rounding marginally 
stable orbit with Keplerian velocity. In the case of a optically 
thin accretion flow like Bondi-type flow in which an angular 
momentum of an accreting matter is smaller than Keplerian 
and radial velocity of the matter is large, it is open to question 
whether a period of a marginally stable orbit determine a min- 
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imum time of variability. Other methods to determine a spin 
parameter of a black hole are proposed. X-ray polarized emis- 
sion from an accretion disk around a black hole (e.g. Conners, 
Stark & Piran 1980), microlensing light curves due to a black 
hole as a lensing object (e.g. Asada, Kasai & Yamamoto 2003) 
and gravitational wave emission from a star orbiting around a 
black hole (e.g. Tanaka et al. 1996, Finn & Thorne 2000) are 
also influenced by a black hole rotation. 

Observation of a black hole shadow is another possible 
method to determine a spin parameter of a rotating black hole 
in future. This is because in the vicinity of a black hole null 
geodesies directly reflect effects of black hole rotation, such 
as frame dragging. As Falcke, Melia & Agol (2000) pointed 
out, in the case of Sgr A* radio interferometers with sufficient 
spatial resolutions may resolve the central black hole shad- 
ows in the optically thin emitting medium by observations of 
sub-millimeter wavelength. A shape of a black hole shadow 
in an optically thin medium are deformed by black hole rota- 
tion (e.g. Bardeen 1973, Chandrasekhar 1983, Falcke, Melia & 
Agol 2000). An apparent shape of such black hole shadow in 
the (jc, y)-plane is described by (Bardeen 1973, Chandrasekhar 
1983) 
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where m and a are a mass of a black hole and an angular mo- 
mentum per unit mass of a black hole, respectively, i is an 
inclination angle between a rotation axis of a black hole and 
a direction to an observer and / is an intervening parameter. 
Here, x and y are normaUzed by a gravitational radius Vg de- 
fined as Vg = Gm/c^, where G is the gravitational constant and 
c is the speed of light. In the case of large inclination ; ^ 90° 
and a large spin parameter a/m ~ 1, a black hole shadow have 
an elongated shape in a direction of a rotation axis of a black 
hole (e.g. Bardeen 1973, Chandrasekhar 1983, Falcke, Melia 
& Agol 2000). If such an elongation feature of a black hole 
shadow in optically thin emitting medium is observed, the spin 
parameter of the black hole can be determined by fitting a con- 
tour of a black hole shadow given by equations (1) and (2) to 
the observed image. 

On the other hand, in the case of a black hole in an opti- 
cally thick and geometrically thin accretion disk, is it possible 
to measure a spin parameter of a black hole from observation of 
a black hole shadow? One of the main purposes of the present 
paper is to answer this. A contour of a black hole shadow 
in an accretion disk is basically calculated by Cunningham & 
Bardeen (1972) and Cunningham (1973, 1975) for a rotating 
black hole. Luminet (1979) and Fukue & Yokoyama (1988) 
calculate apparent images of black hole shadows in accretion 
disks for non-rotating black holes. In the case of a rotating 
black hole, such as a deformed black hole shadow in optically 
thin medium as mentioned above, one might firstly think that 
a spin parameter of a black hole can be determined from a de- 
formed shape of a black hole shadow due to a black hole rota- 
tion. Actually, many authors calculate such deformed shadows 
due to black hole rotation in various contexts (Cunningham & 
Bardeen 1972, Cunningham 1973, 1975, Sikora 1979, Perez & 
Wagoner 1991, Jaroszynski, Wambsganss, & Paczynski 1992, 
Kindl 1995, Quien, Wehrse & Kindl 1996, Begelman & Rees 
1996, Hollywood & Melia 1995, 1997, Bromley, Miller & 
Pariev 1998, Pariev & Bromley 1998, Bromley, Melia & Liu 
2001). In the present paper, we show black hole shadows in ac- 
cretion disks with quite similar shapes for largely different spin 
parameters of central black holes. When we calculate a shape of 
a black hole shadow in an accretion disk, an irmermost radius of 
a luminous part of an accretion disk is significantly important. 
Traditionally, as an emitting disk surrounding a black hole, the 
standard disk (Shakura & Sunyaev 1973, Novikov & Thorne 
1973) is used. In such a disk, the irmermost radius of the lumi- 
nous part of the accretion disk corresponds to a marginally sta- 
ble orbit if we neglect dimming effects due to gravitational red- 
shift in the vicinity of the black hole. Since the surface density 
inside this edge is expected to be extremely smaller than that 
outside the edge, it is usually assumed that there is no luminous 
matter inside the inner edge. However, it has not been proven 
whether the region inside the inner edge is non-luminous or not 
(Fukue 2003). Reynolds and Begelman (1997) discuss in detail 
the possibihty that the material within innermost stable orbit 
may be bright in the X-ray iron line. Recently, from relativistic 
magnetohydrodynamical simulations, Krolik & Hawley (2002) 
proposed that the inner edge of an accretion disk around a black 
hole depends on the definition, i.e. turbulence edge, stress edge, 
reflection edge, and radiation edge are defined by each physical 
process. In the present study, we define the inner edge as an 



innermost radius of a luminous part of an accretion disk. Also, 
in the case of a super-Eddington accretion disk, an optically- 
thick luminous disk extends from a large radial distance to an 
event horizon (e.g. Watarai & Mineshige 2003), i.e. an inner- 
most radius of a super-Eddington accretion disk is smaller than 
a marginally stable orbit. This is because an effective poten- 
tial based on an angular momentum distribution does not have 
a potential minimum for accretion rates, M > 10Le/c^, where 
Le is the Eddington luminosity (Watarai & Mineshige 2003). 

In the present study, we calculate contours of black hole 
shadows in geometrically thin and optically thick accretion 
disks for two kinds of innermost radii, a marginally stable orbit 
and an event horizon. When an innermost radius of an accre- 
tion disk is not close to an event horizon, the radiation from 
the underside of the accretion disks appear within a black hole 
shadow consisting of the direct orbits. Since the intensity of 
such radiation is not strong, we neglect such radiation. In §2, 
calculation methods are shown. We give a qualitative analy- 
sis and a quantitavive analysis of the shape and position of the 
black hole shadow casted by the rotating black hole in the center 
of the accretion disk for two kinds of innermost radii of accre- 
tion disks in §3. We present quite similar black hole shadows 
for largely different spin parameters of black holes in §4. We 
give discussion and conclusions in §5. 

2. CALCULATION METHODS OF BLACK HOLE SHADOWS IN 
ACCRETION DISKS 

In this section, we briefly show a calculation method of a 
contour of a black hole shadow in Kerr spacetime. The propa- 
gation of radiation from a matter around a Kerr black hole is ex- 
tensively studied by Cunningham (1975). The background met- 
ric is written by the Boyer-Lindquist coordinates with c = G = 1 , 
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where A = -2mr + a^, S = r^-i-a^cos^^, and m and a de- 
note the mass and angular momentum per unit mass of the 
black hole, respectively. In the present study, we assume an 
axisymmetric emitting matter distribution with respect to a ro- 
tation axis of a black hole. In such a case, a rotation axis of a 
black hole is equivalent to a rotation axis of an accretion disk. 
In the present study, we only consider a co-rotating accretion 
disk around a rotating black hole. Since we assume the ax- 
isymmetric matter distribution with respect to a rotation axis of 
a black hole, we only need to calculate projected null geodesies 
on the (r, 6) plane in order to calculate a contour of a black hole 
shadow. An integral form of the geodesic equations of (r, 9)- 
componentcan be written as (e.g. Carter 1968, Misner, Thorne 
& Wheeler 1973, Chandrasekhar 1983) 
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where 

eie) = q^ + a^cos^9-p^cot^9, (5) 
R(r) = r^ + {a^-p^-q^)r^ + 2m[q^ + {a-pf]r-q^a^. (6) 
Here, p and q are constants along the geodesic and defined as 
P = LJE, q^ = ClE\ (7) 
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where E, and C are an energy, an angular momentum and 
Carter's constant along the geodesic, respectively. Two impact 
parameters (x, y) in the plane of the sky of an observer can 
be expressed by the two constant, p and q, as (Bardeen 1973, 
Curmingham & Bardeen 1972, Chandrasekhar 1983) 

x = -/7/sin;, (8) 

y = ^^ + cos^;-/5^cot^/, (9) 

where / is an inclination angle between a rotation axis of accre- 
tion disk and a direction to an observer. As several authors have 
done, 9 and r integrals of geodesies can be solved analytically 
by using ElUptic integrals of the first kind and the second kind 
(e.g. Cunningham 1975, Rauch & Blandford 1994). These ana- 
lytic solutions have many different sub-cases depending on the 
values of constants along the trajectory and the spin parameter 
of the black hole. Also, the signs of the equation (4) are cho- 
sen independently (Carter 1968) and depend on the direction 
of the geodesies. Rauch & Blandford (1994) extensively stud- 
ied these analytic solutions and the calculation method of the 
geodesies in a Kerr metric. All patterns of analytic solutions 
of equation (4) are Usted in tables in Appendix A of Rauch & 
Blandford (1994). Note that while Rauch & Blandford (1994) 
give all analytic solutions of r-0 components of the geodesies 
in a Kerr metric, Rauch & Blandford (1994) do not present the 
analytic solutions of the ^-component and the ^-component of 
the geodesies. For details of the calculation methods of the 
geodesies in a Kerr metric, refer to Rauch & Blandford (1994). 
In the present study, we use the analytic formulas of equation 
(4). If we consider an observer at infinity from a black hole, in 
the equation (4) we set r and 9 components of an observer to be 
(ro, ^o) = (oo, 0- When we need to specify ro, we use a value 
ro = 10^ (in units of r^). In the subsequent sections we calcu- 
late black hole shadows in thin accretion disks and we set a 9 
component of an end point of a geodesies to be = 7r/2 in the 
equation (4). In terms of an innermost radius of a luminous part 
of an accretion disk, we consider two types of radii as an exam- 
ple; an event horizon rh+ and a marginally stable orbit r^s. An 
event horizon of a black hole defined as rh+ = m + {np- -a^)^/'^, 
and an expUcit form of a marginally stable orbit r^s can be seen 
in Bardeen, Press & Teukolsky (1972). Scales of black hole 
shadows in the present study are normalized by a gravitational 
radius, Vg, defined as Vg = Gm/c^. 

3. BLACK HOLE SHADOWS IN GEOMETRICALLY THIN 
ACCRETION DISKS 

3.1. Qualitative Features 

In order to understand how a contour of a black hole shadow 
depend on a spin parameter, a/m, of a black hole, an inner edge, 
Tin, of an accretion disk and an inclination angles, we plot 
in figure 1 contours of black holes shadows in optically thick 
accretion disks for a/m = 0, 0.5 and 0.998 (from left panels 
to right panels), r-m = rh+ (thick solid lines) and r^s (thin solid 
lines), and i = 0°, 45° 80° and 89° (from top panels to bottom 
panels). We also plot contours of black hole shadows casted 
by black holes in optically thin emitting matters (e.g. Bardeen 
1973, Chandrasekhar 1983, Falcke, Melia & Agol 2000) by 
dashed Unes in figure 1. These shadows are calculated from 
equations (1) and (2). As a reference frame, in all panels of 
figure 1 , we also plot lines of x = and those of y = by dotted 
lines. The origin, (x, y) = (0, 0), is a position of a black hole, 
i.e. a mass center of a black hole. In this subsection we sum- 
marize quaUtative features of contours of black hole shadows 



in accretion disks. These features are partly investigated in ref- 
erences cited in §1. In subsequent subsections, we investigate 
black hole shadows quantitatively. 

First, we see how a contour of a black hole shadow depends 
on an inclination angle, ;. In the case of ; = 0°, a shape of a 
black hole shadow is circle with its center at the origin, i.e. a 
position of a black hole mass center, for any black hole spin pa- 
rameters because of axisynnmetric luminous matter distribution 
with respect to a rotation axis of a black hole. Although a size 
of a black hole shadow for an optically thin medium plotted by 
a dashed line is nearly independent of a black hole spin parame- 
ter, for both cases of shadows in disks with r^ = rh+ (thick soUd 
lines) and rs„, (thin solid lines) the size of the black hole shadow 
become smaller for a more rapidly rotating black hole. This is 
simply because a radius of an event horizon and a marginally 
stable orbit become smaller for a more rapidly rotating black 
hole. In the case of finite inclination angles, i.e. i = 45°, 80° 
and 89° in figure 1, since accretion disks non-axisymmetrically 
prevent null geodesies from plunging into a black hole horizon, 
black hole shadows for cases = rh+ and deviate from cir- 
cles, and slightly become flat shapes. In the case of a black 
hole shadow for optically thin medium, since there is no occul- 
tation by an accretion disk the shape of the black hole shadow is 
nearly circle even for a finite inchnation angle and a finite value 
of a black hole spin parameter. Especially, in the case of i = 89° 
the black hole shadow in an accretion disk for cases ri„ = rh+ 
or Tsm, underside of a black hole shadow which is seen in the 
case of a shadow in optically thin medium is partly occulted by 
an accretion disk. When an outer radius of an accretion disk is 
not so large, a black hole shadow consists of radiation from a 
backside of an accretion disk can be seen when viewing from 
nearly edge-on, i ~ 90°. In the present study, we assume suffi- 
ciently large accretion disks and do not consider radiation from 
the backside of an accretion disk. 

In the case of = rh+, a horizontal size of a black hole 
shadow becomes larger for a larger inclination angle. In the 
case of a finite inclination angle, in a region of an upper side, 
i.e. y > 0, a length of a geodesies from an observer to a crossing 
point on an accretion disks becomes longer than a geodesies in 
the case of i = 0°. In such cases, some geodesies crossing accre- 
tion disks when ; = 0° directly plunge into black hole horizon, 
and a size of a black hole shadow becomes larger for a larger 
inchnation angle. This effect is strong near a black hole. In the 
case of Tin = a size of a black hole shadow becomes slightly 
larger for a larger inclination angle. But in the case of ri„ = r^s 
this effect is weaker than cases when Hn = rh+. This is because 
when a marginally stable orbit is larger than an event horizon 
for a fixed spin parameter and in the case of r^^ = rb+ geodesies 
consisting contours of black hole shadows are more strongly 
bent than those in the case of = r^s. 

Additionally, in the case of nn = r^s, a nearly maximum in- 
clination angle, i ~ 90°, and a non-rotating black hole, as seen 
in the left bottom panel in figure 1, a black hole shadow has 
two tips near >■ ^ and x ^ ±7. Upper side of a black hole 
shadow is made by geodesies from a back of a black hole and 
these geodesies are strongly bent by a strong gravitational field 
of a black hole. In the region relatively apart from a black hole, 
these bending effects are not so strong and geodesies become 
nearly Newtonian. In such cases, black hole shadows have tips 
which reflect nearly Newtonian geodesies. So, if we put an ac- 
cretion disk with an iimer edge larger than a marginally stable 
orbit and we view such shadows with clear tips when a large 
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inclination angle. 

As mentioned above, in the case of a rotating black hole for 
a larger spin parameter of a black hole an event horizon and 
a marginally stable orbit become smaller and a size of a black 
hole shadows also become smaller. But, this smallness of a 
black hole shadow is not unique properties of effects of black 
hole rotation. For a constant spin parameter of a black hole, 
a size of a black hole shadow becomes smaller for a smaller 
inner edge of an accretion disk. Also, for a constant spin pa- 
rameter of a black hole and a constant inner edge of an accre- 
tion disk, a size of a black hole shadow becomes smaller for a 
smaller black hole mass. What is the unique property of black 
hole rotation? Although in the case of a non-rotating black 
hole a black hole shadow is perfectly symmetric with respect 
to a rotation axis of an accretion disk, as shown in figure 1, in 
the case of a rotating black hole a black hole shadow becomes 
non-axisymmetric with respect to a rotation axis. This is be- 
cause geodesies near a black hole suffers from effects of frame 
dragging and such a geodesic becomes non-axisymmetric with 
respect to a rotation axis. This non-axisymmetric shape of 
the black hole shadow can be also seen in past studies (Cun- 
ningham & Bardeen 1972, Cunningham 1973, 1975, Sikora 
1979, Perez & Wagoner 1991, Jaroszynski, Wambsganss, & 
Paczyhski 1992, Kindl 1995, Quien, Wehrse & Kindl 1996, 
Begelman & Rees 1996, Hollywood & Melia 1997, Bromley, 
Miller & Pariev 1998, Pariev & Bromley 1998, Bromley, Melia 
& Liu 2001). If luminous matters around a black hole are dis- 
tributed non-axisymmetrically, the non-axisymmetric property 
of a black hole shadow is not unique effect of black hole rota- 
tion. 

3.2. Maximum Size of a Black Hole Shadow 

In this section, we see quantitative features of the maximum 
size of a black hole shadow. In figure 2, we show dependence of 
the maximum size of a black hole shadow on a spin parameter 
of a black hole and an inclination angle for three kinds of sur- 
rounding mediums; an accretion disk with innermost radius of 
Hn = rb+ (dotted lines) and an accretion disk with = r^s (solid 
lines), and optically thin emitting medium (dash-dotted lines). 
Note that the maximum size of a black hole shadow generally 
does not coincide with the maximum horizontal size of a black 
hole shadow because of black hole rotation. In figure 1, this 
feature can be clearly seen in black hole shadows in disks for 
a/m = 0.998andj = 80°. 

We can see gravitational focusing effects of geodesies in fig- 
ure 2. For example, in the case of a non-rotating black hole 
a radius of a marginally stable orbit is 6. So, a diameter of a 
circle of a marginally stable orbit is 12. However, for an incli- 
nation angle i = 0°, ri„ = and a non-rotating black hole the 
maximum size of a black hole shadow is nealy 13.85 in figure 
2 because of gravitational focusing effects. 

For any spin parameter of a black hole, the maximum size 
of a black hole shadow becomes larger for a larger inclination 
angle. The reason for this is explained in the previous section. 
For any spin parameter of a black hole and a fixed inclination 
angle a black hole shadow in an accretion disk with = rt+ is 
smallest among black hole shadows for these three kinds of sur- 
rounding medium. In the case of a slowly rotating black hole 
or a non-rotating black hole (a/m < 0.5), the maximum size 
of a black hole shadow around optically thin medium is among 
those in the cases of the other two disks for same inclination 
angle. On the other hand, in the case of a rapidly rotating black 



hole (a/m > 0.8) the maximum size of a black hole shadow 
around optically thin medium is larger than the black hole shad- 
ows in other two kinds of accretion disks for a fixed incUnation 
angle. Although the maximum size of a black hole shadow in 
optically thin medium concentrates on a narrow range of a size 
(i.e. from ~ 9.6 to ~ 10.4), the maximum size of a black hole 
shadow in an accretion disk expands a relatively wide range of 
a size. This is because the event horizon and marginally stable 
orbit depend on a spin parameter of a black hole as mentioned 
above. 

In the case of an accretion disk with = rh+, the maximum 
size of a black hole shadow is ~ 8.7 for i = 90° and a/m = 0. 
On the other hand, the minimum size of a black hole shadow 
is ~ 4.0 for / = 0° and a/m= 1. When an accretion disk with 
''in = '"ms, the maximum size of a black hole shadow is ^ 14.7 
for i = 90° and a/m = 0, and the minimum size of a black hole 
shadow is ~ 4.2 for J = 0° and a/m = 1. In contrast to the case of 
a black hole shadow in optically thin medium for which a size 
of a black hole shadow is in a narrow range as mentioned above, 
if an inner edge of an accretion disk is equal to or smaller than 
a marginally stable orbit, a size of a black hole shadow expands 
in a wider range from 4.0 to 14.7 than a black hole shadow in 
optically thin medium as shown in figure 2. This is a remark- 
able quantitative feature of a black hole shadow in an accretion 
disk. 

3.3. Flatness of Black Hole Shadow 

Here, for a single black hole shadow we consider a ratio of 
the minimum width of a black hole shadow to the maximum 
width of a black hole shadow. This ratio is important because 
this ratio is a good indicator of an inclination angle between 
a rotation axis of a black hole and a direction of an observer. 
This ratio is similar to a flatness of a black hole shadow. This 
ratio is a dimensionless quantity. Note that a line segment of 
the minimum width and a line segment of the maximum width 
generally do not cross perpendicularly. In figure 3, we show 
dependence of the ratio on a spin parameter of a black hole for 
three inclination angles (45°, 80° and 89°) and two kinds of 
accretion disks with = rh+ (dotted lines) and rin = r^^s (solid 
lines). In the case of an inchnation angle of i = 0°, this ratio 
is unity. Although the ratio strongly depends on an incUnation 
angle, the ratio weakly depends on a spin parameter of a black 
hole. In the both cases of black hole shadows in accretion disks 
with Tin = rh+ and ri„ = rms, for an inclination angle / = 89° this 
ratio is nearly 0.5 for most of a range of a spin parameter. This 
value of ~ 0.5 is a lower limit of the ratio. 

3.4. Horizontal Shift of a Black Hole Shadow from a Rotation 
Axis of a Black Hole 

As mentioned in §3.1, if we consider an axisymmetric lu- 
minous matter distribution with respect to a rotation axis of a 
black hole, a horizontal shift of a black hole shadow, which is 
a displacement of a shadow center in the direction of x axis, 
shows a significant unique property of effects of a black hole 
rotation. In order to quantify this shift, we define x coordinates, 
Xc, asxc = (xniax+.^min)/2 where Xmax and Xmin are the maximum 
x-coordinate and the minimum x-coordinate of a contour of a 
black hole shadow, respectively. In the case of a non-rotating 
black hole, a shape and a position of a black hole shadow are 
synamettic with respect to a rotation axis of a black hole, i.e. 
y axis. So, in this case, Xc = 0. On the other hand, in the case 
of a rotating black hole a shape and a position of a black hole 
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shadow is not symmetric with respect to a rotation axis of a 
black hole. So, in this case a value of Xc is finite. This is be- 
cause frame dragging effects of geodesies due to a black hole 
rotation break the symmetric property of a black hole shadow 
with respect to a rotation axis. 

In figure 4, we show the shift of the center of the black 
hole shadow in the direction of x axis for three inclination an- 
gles (45°, 80° and 89°) and two kinds of accretion disks with 
rin = i~h+ (dotted lines) and rin = r^s (solid lines). In this fig- 
ure, we can see that an extent of the horizontal shift is roughly 
proportional to a spin parameter of a black hole for a fixed in- 
cUnation angle. An extent of the maximum horizontal shift is 
about 1 .45 in the case of the maximum rotation, the maximum 
inclination angle and rin = ?"h+. A horizontal shift for the case 
of Tin = rh+ is larger than that for the case of ri„ = r^s for same 
inclination angle and same spin parameter of a black hole. This 
is because effects of frame dragging due to black hole rotation 
are more significant in a nearer region of a black hole. These 
strong frame dragging effects result in a large shift of a black 
hole shadow in x direction. If we can measure this horizontal 
shift observationally, we can well constrain a spin parameter 
of a black hole. For example, in the case of an accretion disk 
with an inner edge between an event horizon and a marginally 
stable orbit, a spin parameter of a black hole is constrained be- 
tween 0.65 and 0.83 for an incUnation angle of i = 45° and a 
horizontal shift of 0.3. Here, we assume a geometrically thin 
accretion disk. In order to measure a horizontal shift of a black 
hole shadow, which is a good indicator of a spin parameter of 
a black hole, we need to know a position of a rotation axis of 
a black hole and a position of a center of a black hole shadow. 
However, it seems to be difficult to determine a position of a 
rotation axis of black holes observationally. So, although a hor- 
izontal shift of a black hole shadow is a good indicator of a spin 
parameter of a black hole, this method is not a practical method 
to determine a spin parameter of a black hole by observation 
of a black hole shadow. If we can determine a position of a 
rotation axis of a black hole by some methods, this method is 
practical for determination of a spin parameter of a black hole. 
In the next section, we consider another method to determine a 
spin parameter of a black hole from an image of a black hole 
shadow. 

4. CAN WE MEASURE A SPIN PARAMETER OF A BLACK 
HOLE IN AN ACCRETION DISK BY OBSERVATION OF A 
BLACK HOLE SHADOW? 

4.1. Quite Similar Black Hole Shadows for Largely Different 
Spin Parameters of Black Holes 

In principle, we can determine a spin parameter of a black 
hole a /m, an inclination angle i and an innermost radius of an 
accretion disk by fitting a calculated contour of a black hole 
shadow to an observed image of a black hole shadow obtained 
by observational facilities with sufficient spatial resolutions and 
sufficient detection Umits even in the case that the position of 
a black hole shadow are not known. However, it is practically 
difficult to determine spin parameters of black hole shadows by 
this fitting method. This is because there are black hole shad- 
ows with a quite similar shape and size for largely different 
spin parameters of black holes. In figure 5, we show examples 
of such black hole shadows with nearly same maximum width 
of a black hole shadow and nearly same ratio of the minimum 
width to the maximum width of a black hole shadow. The max- 
imum widths of the black hole shadows are ~ 6.2 for black 



hole shadows in panels of (a), (b) and (c), and ~ 7.7 for black 
hole shadows in panels of (d), (e) and (f), respectively. The 
flatness ratios are ~ 0.86 for black hole shadows in panels of 
(a), (b) and (c), and ~ 0.60 for black hole shadows in panels of 

(d) , (e) and (f), respectively. Here, the adapted parameters are 
(a) a/m = 0, i = 45° and = rh+, (b) a/m = 0.5, i = 45° and 
rin = 1.05rh+, (c) a/m = 0.958,?' = 45° andri„ = rn,s, (d)a/m = 0, 
i = 80° and = rh+, (e) a/m = 0.5, i = 80° and r^ = 1.05rh+ 
and (f) a/m = 0.955, i = 80° and rin = r^s. In each panel of 
figure 5, a position of a black hole mass center is plotted by a 
white dot. a rotation axis of an accretion disk is plotted by a 
gray line. In the case of a rotating black hole, this line is also 
a rotation axis of a black hole. A horizontal direction in each 
panel of figure 5 represents a direction of the maximum width 
of a black hole shadow. As shown in figure 1 and figure 4, in 
the case of a finite spin parameter of a black hole, a black hole 
shadow is non-axisymmetric with respect to a rotation axis of 
a black hole. In figure 5, extents of shifts of centers of black 
holes from rotation axes of black holes are (b) 0.22, (c) 0.38, 

(e) 0.50 and (f) 0.95, respectively. As mentioned above these 
shifts are nearly proportional to spin parameters of black holes 
(see figure 4). In the case of a rotating black hole a direction 
of a segment of the maximum width of a black holes is not per- 
pendicular to a rotation axis of a black hole. Thus, as panels 
of (b), (c), (e) and (f) show, rotation axes of black holes are in- 
cUned from vertical directions. These inclination angles from 
vertical directions are (b) 6.2°, (c) 12°, (e) 4.2° and (f) 9.2°, 
respectively. One of the largely different points between black 
hole shadows with a similar shape and a size is a distance be- 
tween a rotation axis of an accretion disk and a center of a black 
hole shadow. If we know a position of a rotation axis of a black 
hole, a spin parameter of a black hole can be well constrained 
from an image of a black hole shadow. But, as mentioned in the 
previous section, it seems difficult to determine a position of a 
rotation axis of a black hole. 

Here, we newly introduce a bisector axis of a black hole 
shadow named a shadow axis. We define this shadow axis as a 
bisector perpendicular to a line segment of the maximum width 
of a black hole shadow. In each panel of figure 5, we plot this 
shadow axis by a white line in the case of a rotating black hole. 
In the case of a non-rotating black hole, a shadow axis coincides 
with a rotation axis of a black hole shadow. We can easily de- 
termine a position of the shadow axis by observation of a black 
hole shadow in an accretion disk. In figure 5, the minimum 
intervals between a mass center of a black hole and a shadow 
axis are (b) 0.32, (c) 0.58, (e) 0.59 and (f) 1.1, respectively. 
This minimum interval strongly depends on a spin parameter of 
a black hole and can be measured if a position of a mass center 
of a black hole is known and a black hole shadow is observed. 
So, in the next section, we investigate a relation between a spin 
parameter of a black hole and the minimum interval between a 
mass center of a black hole and a shadow axis. 

4.2. Minimum Interval between a Black Hole Mass Center 

and a Shadow Axis 

In figure 6, we plot the minimum interval between a black 
hole mass center and a shadow axis for all range of a spin pa- 
rameter of a black hole. Figure 6 is similar to figure 4, but there 
are differences between lines in figure 4 and those in figure 6. 
In figure 6, we plot the minimum intervals for three inclination 
angles (45°, 80° and 89°) and two kinds of accretion disks with 
r-m = rb+ {dotted lines) and ria = r^s {solid lines). As mentioned 



6 



Rohta Takahashi 



in §3.1, for a non-rotating black hole around an axisymmet- 
ric luminous accretion disk with respect to a rotation axis of a 
black hole, a black hole shadow is symmetric with respect to a 
rotation axis of a black hole. Thus, in figure 6 for a non-rotating 
black hole, i.e. a/m = 0, the minimum interval is null. On the 
other hand, in the case of a rotating black hole, i.e. a/m> 0, 
the shape and position of the black hole shadow is not symmet- 
ric with respect to a rotation axis of a black hole. So, in such 
case the minimum interval is finite. In figure 6, we can see that 
an extent of the minimum interval is roughly proportional to a 
spin parameter of a black hole for a fixed inchnation angle. The 
maximum extent of the minimum interval is about 1 .5 in the 
case of the maximum rotation, the maximum inclination angle 
and Tin = rh+. The minimum interval for the case of = rh+ 
is larger than the minimum interval for the case of = r^s 
for same inclination angle and same spin parameter of a black 
hole. This is because effects of frame dragging due to black 
hole rotation are more significant in a nearer region of a black 
hole, as mentioned in the previous section. If we can measure 
such minimum interval, a spin parameter of a black hole can 
be well constrained. But if not, quite subtle fitting of a contour 
of a black hole shadow to an observed image of a black hole 
shadow is required in order to determine a spin parameter of a 
black hole and other physical parameters from an observed im- 
age of a black hole shadow. For example, for an accretion disk 
with an inner edge between an event horizon and a marginally 
stable orbit with an inclination angle of ! = 45° and a horizontal 
shift of a black hole shadow to be 0.5, from figure 6 a spin pa- 
rameter of a black hole is constrained between 0.75 and 0.89. 
Thus, in order to measure a spin parameter of a black hole it is 
crucially important to determine a position of black hole mass 
center in a region of a black hole shadow. Required positional 
accuracy to determine a spin parameter of a black hole is lower 
than about 1 .5 (in units of rg). This is the case of the maximum 
shift of black hole shadow center from a rotation axis. In the 
last section, we discuss determination methods of a mass center 
of a black hole. 

Note that in order to determine a spin parameter of a black 
hole from an image of a black hole shadow we do not always 
need to determine a position of a black hole mass center. For 
some cases, we can determine a spin parameter of a black hole 
by fitting a calculated contour of a black hole shadow to an 
observed image of a black hole shadow. One of the examples 
of such cases is a black hole shadow with tips which can be 
seen when a relatively large inner edge of an accretion disk as 
mentioned in §3.1. In bottom panels of figure 1, in the case 
of Tin = there are tips in contours of black hole shadows 
for a/m = and 0.5. When a/m = there are two tips, when 
a/m = 0.5 there is a single tip and when a/m = 0.998 there is no 
tips. So, if there are one or two tips in contours of black hole 
shadows, fitting of a calculated contour of a black hole shadow 
to an observed black hole shadow is more easy than the cases 
with no tips in a contour of a black hole shadow. But, the black 
hole shadow with tips is a quite limited case. 

5. DISCUSSION AND CONCLUSIONS 

5.1. How is the Position of the Center of Mass of a Black Hole 

Determined? 

As mentioned in the previous section, in order to measure 
a spin parameter of a black hole from an observed image of a 
black hole shadow it is important to determine a position of a 
black hole mass center. In the case of Sgr A*, stellar motions 



around black holes are actually observed (e.g. Genzel et al. 
2003). From these stellar motions, in principle we can deter- 
mine the position of a mass center of a black hole. Astrometry 
satellite missions with positional accuracy of 10 /x arc second 
which will give a plenty of data of stellar motions around Sgr 
A* are planned: GAIA for optical observation (e.g. Belokurov 
& Evans 2002 and references therein) and JASMINE (Japan 
Astrometry Satellite Mission for INfrared Exploration) (Gouda 
2002) for infrared observation in which effects of interstellar 
absorption and scattering are not very significant compared to 
those of optical observation. 

There is a simple additional method to determine a black hole 
mass center. Distribution of brightness of the radiation from the 
accreting material around the black hole are influenced by the 
Doppler and gravitational boosting effects. For example, Lu- 
minet (1979) provides detailed derivations and discussions of 
the brightness distribution over the image of the accretion disk. 
On the other hand, an angular momentum, u), of the frame drag- 
ging due to a black hole rotation is described as. 



u =■ 



St<t> 



2mra 



(10) 



(11) 



(r^ + a^)^-a^Asm^e' 
In the equatorial plane, this angular momentum is expanded as 
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Because the frame dragging effect decreases in proportion to 
ajr^ as shown here for a sufficiently large radial distance, r, 
from a mass center of a black hole, and the gravitational po- 
tential decreases only proportional to 1 jr, the frame dragging 
effects on the propagation of photons become negligible at a 
sufficient distance from the black hole. The image at these large 
radii is formed only by the influence of Doppler effect from the 
moving gas and gravitational redshift in a spherically symmet- 
ric potential. This image should be clearly visible, since the 
shadow of the black hole is assumed to be clearly visible. If the 
motions and density distribution of the gas are axisymmetric, 
than the image at large radii would be the same as produced 
by a non-rotating black hole. One can then use calculations by 
Luminet (1979) to fit the observed brightness distribution and 
determine the center of mass and the position of the rotation 
axis of the black hole as the parameters of this fit. Then, it 
makes sense to use both figure 4 and figure 6 to constrain the 
spin parameter of the black hole. 

5.2. Observational Feasibility of Black Hole Shadows 

Although in the present study we investigate black hole shad- 
ows in accretion disks from mainly academic interests as sev- 
eral past studies about black hole shadows in accretion disks 
(e.g. Luminet 1979, Fukue & Yokoyama 1988, Fukue 2003), 
observational feasibiUties of black hole shadows must be con- 
sidered. Falcke, Agol & Melia (2000) shows that a black hole 
shadow in Sgr A* can be observed at sub-millimeter wave- 
lengths at which a scattering radius due to interstellar medium 
is sufficiently small in order to observe a black hole shadow. 
Also, in nearby elliptical galaxy M87 there is a slightly smaller 
black hole shadow than that in Sgr A*. Both of observed spec- 
trum of SgrA* and M87 will be well reproduced by advection 
dominated accretion flows (ADAFs) or recently radiatively in- 
efficient accretion flows (RIAFs) for Sgr A*. In such accretion 
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models, matter densities are quite low and it is expected that 
optically thick accretion disks are not formed. So, black hole 
shadows in Sgr A* and M87 seems to be seen as black hole 
shadows not in accretion disks but in optically thin medium as 
Falcke, Melia & Agol (2000) have calculated. In the case of 
M87, Junor, Viretta & Livio (1999) report observations at 43 
GHz of inner regions of M87 and shows a remarkably broad jet 
near the center and strong colUmation of the jets occurring at 
^30-100 Schwarzschild radii from a black hole. These results 
are consistent with the hypothesis that jets are formed by an ac- 
cretion disk around the central black hole. In the case of M87 
we can not completely deny the observational feasibiUty of a 
black hole shadows in an accretion disk by next-generation ra- 
dio interferometers. In such cases, the present study about black 
hole shadows in accretion disks is applicable to the black hole 
shadow in M87. The central region of M87 is one of the ob- 
servational targets of the planned next-generation VLBI Space 
Observatory Program (VSOP-2) (Hirabayashi et al. 2001) with 
extremely high resolution. Theoretically, the black hole shadow 
in the center of M87 need more studies which include radiative 
transfer calculations with Synchrotron source functions of ac- 
cretion disks and jets. 

5.3. Conclusions 

In the present study, we make a qualitative analysis and a 
quantitative analysis of the shape and position of the black hole 
shadow casted by a rotating black hole on the optically thick 
accretion disk and its dependence on the angular momentum 
of the black hole. We consider two types of innermost radii, 
an event horizon and a marginally stable orbit, of an accretion 
disk. There are black hole shadows with quite similar sizes and 
shapes for largely different spin parameters of black holes and 
same masses of black holes. In the present study, we newly in- 
troduce a bisector axis of a black hole shadow named a shadow 



axis. We define this shadow axis as a bisector perpendicular to 
a line segment of a maximum width of a black hole shadow. 
We can determine a position of shadow axis of a black hole 
shadow by observation of a black hole shadow. The minimum 
interval between a black hole mass center and a shadow axis 
is important. This is because the minimum interval is quite 
roughly proportional to a spin parameter of a black hole for a 
fixed incUnation angle. The maximum extent of the minimum 
interval is about 1 .5 in the case of the maximum rotation, the 
maximum inclination angle and = rh+. If we can measure 
such minimum interval, i.e. if we can determine positions of a 
black hole mass center and a shadow axis, a spin parameter of 
a black hole can be well constrained. Determination methods 
of a mass center of a black hole are discussed in the previous 
subsection. 

Lastly, we adress the assumptions used in the present study. 
In the present study we only consider a geometrically thin ac- 
cretion disk and calculate direct geodesies when calculating 
black hole shadows. An accretion disk with finite thickness 
may cause self-eclipsing effects of an accretion disk for large 
inclination angle and a black hole shadow may be deformed. 
So, this self-eclipse effect could be important for high incUna- 
tion angles and need study in future. There are also geodesies 
which come from an underside of accretion disks. These pho- 
tons also changes shapes of black hole shadows. This effect 
also need study in future. 

The author is grateful to the referee, V. I. Pariev, for helpful 
comments on discussion, K. Ohsuga, K. Watarai and Y. Kato 
for their useful discussions, H. Hirabayashi, H. Kamaya and S. 
Mineshige for their useful comments and continuous encour- 
agement. This work was supported in part by the Grant-in Aid 
of the Ministry of Education, Culture, Sports, and Science and 
Technology. 
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Fig. 1 . — Examples of contours of black hole shadows in geometrically thin optically thick accretion disks with innermost radii, Hn, of = ri,+ {thick solid lines) 
and rju = rms (thin solid lines) and in optically thin emitting medium (dashed lines). Lines of x = and >> = are plotted by dotted hnes. A projected position of a 
black hole is (x, y) = (0, 0). Projected lines of rotation axes of accretion disks or those of black holes are x = 0. In the cases of a finite spin parameter of a black hole 
and a finite inclination angle, we can see an axisymmetric shape and position of a black hole shadow with respect to a rotation axis of an accretion disk. 
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Fig. 2. — Relations between a spin parameter of a black hole and the maximum size of a black hole shadow for four values of inclination angles (0°, 45°, 80° and 
89°) and three kinds of surrounding mediums; accretion disks with innermost radius of = (dotted lines) and those with rj^ = rms (solid lines), and in optically 
thin emitting medium (dash-dotted lines). 
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Fig. 3. — Ratios of the minimum size of a black hole shadow to the maximum size of a black hole shadow for three inclination angles (45°, 80° and 89°) and 
two kinds of accretion disks with rjn = rt^+ (dotted lines) and rj], = rms (solid lines). In the case of an inclination angle of / = 0°, this ratio is unity. The ratio strongly 
depends on an inclination angle between a rotation axis and a direction of an observer, and weakly depends on a spin parameter of a black hole. 
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Fig. 4. — Shifts of centers of black hole shadows in the direction of x axes from rotation axes for three inclination angles (45°, 80° and 89°) and two kinds of 
accretion disks with ri„ = ri,+ (dotted lines) and ri„ = rn,s (solid lines). As this figure shows, the shift is nearly proportional to a spin parameter of a black hole. In 
the case of an inclination angle of i = 0° , an extent of this shift is null. If we can know an extent of a horizontal shift of a black hole shadow from a rotation axis of 
a black hole and an inclination angle, a spin parameter of a black hole is well constrained. 
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Fig. 6. — Relations between a spin parameter of a black hole and the minimum interval between a mass center of a black hole and a shadow axis for three 
inclination angles (45°, 80° and 89°) and two kinds of accretion disks with Tin = rh+ (dotted lines) and rin = rms {solid lines). 



